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ABSTRACT 


Presented  in  the  report  are:  some  results  of  machine  computa¬ 
tions  of  the  scattering  functions,  the  extinction-,  scattering-,  and 
absorption-coefficients  for  pure  absorbing  carbon  particles,  and  the 
compound  carbon  particles  with  concentric  water  shell  of  varying 
thickness;  and  an  analysis  of  the  absorption  properties  of  various 
single  particles. 

The  scattering  and  absorption  properties  of  single  particles  are 
used  to  calculate  the  absorptivity,  transmissivity,  and  reflectivity 
of  three  cloud  models:  a  pure  carbon  particle  cloud,  a  cloud  of  water 
droplets  and  carbon  particles,  and  a  cloud  of  compound  concentric  carbon- 
water  particles.  All  three  clouds  are  1,000  meters  thick  and  near  the 
ground. 

The  temperature  change  of  one  particle  caused  by  radiation  absorp¬ 
tion  has  been  computed,  and  the  total  energy  for  evaporation  of  a  cloud 
mass  has  been  estimated. 

From  these  calculations  it  can  be  concluded  that  a  total  carbon 
mass  of  the  order  of  several  hundred  pounds  per  kn»  will  be  necessary 
to  evaporate  a  stratus-type  cloud. 
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THEORETICAL  CONSIDERATIONS 
ON  THE  EFFECTIVENESS  OF  CARBON  SEEDING 


INTRODUCTION 

The  first  big  step  in  artificial  cloud  modification  was  taken  in  19k6 
when,  through  experiments  conducted  by  Schaefer-*-  and  Vonnegut,  the  artifi¬ 
cial  stimulation  of  precipitation  from  supercooled  clouds  became  feasible. 
The  teghnique  applied  by  the  two  experimentalists  was  based  on  FLndeisen’s 
theory0  that  in  supercooled  clouds  the  Introduction  of  nuclei  could  initiate 
the  formation  of  precipitation  particles.  Schaefer1-  found  that  small  pieces 
of  dry  ice  could  form  a  large  number  of  small  ice  crystals  in  such  a  cloud, 
and  Vonnegut  discovered  that,  also,  tiny  crystals  of  silver  iodide  would  be 
very  active  as  sublimation  nuclei.  Many  successful  seeding  experiments  with 
supercooled  cumulus  clouds  as  well  as  supercooled  stratus  clouds  have  since 
been  carried  out. 

However,  both  of  the  above-mentioned  seeding  methods  are  confined  to 
the  application  for  supercooled  clouds  only,  whose  temperature  is  below  0°C; 
and  all  attempts  to  find  a  reliable  method  for  the  dissipation  of  warm 
clouds  and  fogs  have  been  unsuccessful.  The  principal  difficulty  with  warm- 
cloud  seeding  is  that  there  is  no  "latent  energy"  in  the  cloud  available 
which,  after  liberation,  could  carry  on  the  dissipation  process  by  itself. 

In  the  case  of  supercooled  clouds,  this  latent  energy  is  the  energy  of 
fusion  which  becomes  free  after  seeding  and  which  has  been  provided  in  an 
indirect  way  by  the  decrease  of  the  entropy  of  the  system  by  nature  itself. 
Whatever  energy  for  a  precipitation  process  or  evaporation  process  is  needed 
in  a  warm  cloud,  therefore,  has  to  be  fully  supplied  from  an  outside  source. 
As  such  an  energy  source,  the  radiation  energy  from  the  sun  is,  of  course, 
the  most  apparent.  There  ar.e  other  forms  of  energy  which  could  possibly  be 
used;  for  instance,  electrostatic  forces  for  dissipation  of  fogs  and  clouds. 
These,  however,  will  not  be  discussed  in  this  report. 

Using  the  properties  of  carbon-black  particles  suspended  in  the  atmos¬ 
phere,  several  attempts  have  been  made  to  absorb  and  emit  sunlight,  thus 
influencing,  artificially,  the  energy  processes  in  the  atmosphere  and 
especially  in  clouds.  In  recent  years,  several  investigators  have  made 
estimates  of  the  energy,  absorbed  in  carbon  spheres  and  clouds  of  carbon  par¬ 
ticles;  and  various  hypotheses  have  been  developed  on  the  processes  which 
might  occur  in  carbon-seeded  clouds.  However,  all  these  computations  of  the 
energy  absorbed  by  carbon  particles  have  been  based  on  the  so-called  bulk 
method  which  uses  the  equivalent  solid-mass  layer  for  its  absorption  calcu¬ 
lations,  or  they  were  based  on  very  rough  estimates  of  the  absorption  cross 
section.  Such  results,  of  course,  may  be  ponsiderably  different  from  exact 
solutions  of  the  formulas  of  the  Mie  theory. ^ 


Actual  carbon- seeding  experiments  have  been  carried  out  by  Van  Straten, 
Ruskin,  Dinger,  and  Masteribrook^  and  by  personnel  of  the  Air  Force  Geophysics 
Research  Directorate,  but  none  of  these  experimental  tests  of  carbon  seeding 
gave  a  decisive  answer  to  the  question  of  the  effectiveness  of  this  method. 

It  therefore  seemed  desirable  to  conduct  some  theoretical  studies  on  the  feasi¬ 
bility  of  carbon  seeding,  including  an  evaluation  of  the  absorption  properties 
of  carbon- cloud  elements  according  to  the  Mie  theory**  and  to  find  from  this  a 
sound  basis  for  future  experiments. 

DISCUSSION 

Absorption  and  Scattering  Properties  of  Single  Carbon  and  Carbon-Water 
Particles 


As  mentioned  previously,  for  a  comprehensive  solution  of  the  problem  it 
is  necessary  to  investigate  the  optical  properties  of  pure  carbon  particles 
as  well  as  compound  carbon  and  water  particles.  The  compound  carbon-water 
particle,  a  carbon  sphere  with  a  concentric  water  shell  around  it,  will  be  an 
idealization  to  some  extent  and  will  probably  occur,  in  this  form,  only  in 
the  case  of  condensation  of  water  around  the  carbon  center.  However,  it  is 
the  only  combination  carbon-water  droplet  which  can  be  treated  theoretically. 
The  exact  formulas  for  such  compound  particles  have  been  worked  out  by  Aden 
and  Kerker°  and  are  an  extension  of  Mie's  theory.*1  They  consider  three  dif¬ 
ferent  media  separated  by  two  boundary  layers.  With  this  theory,  the  absorp¬ 
tion  and  scatterxng  properties  of  pure  spherical  homogeneous  particles  with 
complex  refractive  index,  and  also  of  spherical  particles  with  concentric 
shell  of  varying  thickness  can  be  computed.  The  following  notations  will.be 
used: 


9  =  Scattering  angle  (9  =  0  is  forward  scattering) 

a  =  Radius  of  the  inner  sphere  (carbon  particle) 

b  =  Outer  -  radius  of  the  shell 

=  Wavelength  of  the  incident  light 
m  =  n  (l  -  n'i)  =  complex  refractive  index 
ot  =  2  Tf  a/  =  size  parameter  for  a 

'f'  =  2  IT  la/  /\  =  size  parameter  for  b 

b/a  =  Specific  ratio. 

According  to  Mie's  theory  and  from  Aden  and  Kerker's  treatment,^  the  fol¬ 
lowing  quantities  are  derived. 


The  cross  sections  for  extinction  and  scattering, 
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The  functions  Sq  and  bQ  depend  on  the  refractive  index  and  the  size  parame¬ 
ter^);  the  lTn  and  f  n  are  functions  of  the  scattering  angle  only. 

Numerical  computations  have  been  carried  out  for  Qg,  Qs,  i^,  and  ig  for 
pure  carbon  particles,  refractive  index  m  =  1.59  -  0.66 I  (b/a  =  l),  and  car¬ 
bon  spheres  with  a  concentric  water  shell  (b/a  variable)  for  size  parameter 
values  from  0.1  to  200.  The  computations  have  been  carried  out  under 
Signal  Corps  contract  on  the  IBM  704  computer  of  the  National  Bureau  of  Stand¬ 
ards,  Washington,  D.  C. 

In  Figs.  1  through  7  the  computed  total  extinction-,  scattering-,  and 
absorption-coefficients  for  compound  water-carbon  particles  for  various 
ratios  b/a  from  1  to  0  are  plotted  as  functions  of  the  size  parameter.  Since 
per  definition,  the  Qg,  Qg,  and  Qa  give  the  ratio  between  an  area  with  which 
a  particle  interferes  with  the  incident  light  and  the  actual  geometrical 
cross  section,  it  can  be  seen  that  the  effective  absorption  cross  section  of 
a  pure  carbon  particle  (Fig.‘  1)  for  1  <  a  <  4  becomes  almost  1  1/2  its  geomet¬ 
rical  cross  section.  For  a  4.  0.5  the  extinction  is  practically  due  to 
absorptic  only,  and  for  a  >  2  scattering  and  absorption  coefficients  are 
of  the  same  order  of  magnitude. 

Tbr  absorbing  particles,  the  Q- functions  are  very  smooth  and  do  not  show 
the  oscillations  which  are  so  characteristic  for  the  extinction  coefficient 
of  nonabsorbing  particles.  However,  as  the  thickness  of  the  water  shell 
around  the  carbon  particle  increases,  more  and  more  of  the  effect  of  the 
damping  of  the  imaginary  part  of  the  refractive  index  is  lost,  and  the  oscil¬ 
lations  in  the  Q- functions  become  stronger  (b/a  -  1.1,  1.2,  1.5,  2,  5  and  10 
in  Figs.  2  through  7,  respectively)  until  for  very  large  b/a  ratios  the  Q- 
functions  approach  the  values  for  pure  water  droplets.* 


*The  curves  in  Figs.  5  through  7  for  values  ^  }  10  do  not  contain  enough 
computed  points  to  identify  them  unambiguously.  Qg  ( y)  for  b/a  -  10  is 
practically  already  equal  to  Qs  (oc)  for  pure  water. 
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Extinction-,  scattering-,  ana  absorption-coefficients 
for  carbon  particles  of  varying  water  shell  thickness. 
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In  Jig.  8,  the  absorption  coefficient  ^  (referring  to  a  pure  carbon 
particle  as  100  percent),  has  been  plotted  as  a  function  of  b/a  for  constant 
total  particle  size  (  f'  =  const).  It  can  be  seen  from  this  graph  that  for 
larger  particles  ('&'  >10)  the  absorption  coefficient  of  a  carbon  particle 
will  be  increased  even  by  a  thin  water  shell  (b/a  =  1.1).  However,  if  b/a 
increases  further,  the  water  shell  begins  to  reduce  the  absorption  effect 
of  the  carbon  nucleus.  This  shielding  effect  of  the  water  shell  is  the 
stronger,  the  smaller  the  total  particle  size  parameter. 

According  to  the  definition  for  C^,  the  absorption  cross  section  of 
one  particle  is 


ca  =  r  •  7T  .  Qa, 


(7) 


where  r  is  the  radius  of  the  outer  shell;  i.e.,  in  the  case  of  compound 
particles,  b  is  used  instead  of  r.  Prom  this  it  follows  that  the  absorption 
effect  Ca  of  a  compound  particle  which  grows  bigger  and  bigger  disappears 
only  if  Qa(b)  decreases  more  than  b^  increases  with  increasing  b.  However, 
since  Qg  is  a  function  of  oc  and  ,  or  of  either  «  or  ^  plus  b/a,  differ¬ 
ent  results  are  obtained  for  equation  (7),  depending  on  whether  b/a  or  ^ 
is  kept  constant.  If  and  P'  grow  proportionally  (b/a  =  const),  the  increase 
of  b2  will  always  be  of  higher  order  than  the  one  for  Qg,  and  therefore  Ca 
— *  00  as  ^  increases  (Fig.  9)*  'Whereas,  if  oc  i8  kept  constant  and  the 
water  is  allowed  to  grow  around  the  carbon  particle  (i.e.,  b/a  variable),  CL 
(  *)  will  decrease  much  more  for  increasing  ^(Figs.  10  and  11).  In  fig.  8 
it  is  seen  that  even  the  efficiency  factor  for  absorption  C^  itself  increases 
above  the  one  for  a  pure  carbon  particle  if  a  thin  water  shell  forms  around 
the  caybon  core.  In  figs.  10  and  11  for  cC  =0.1  and  ct  =  1,  the  function 
Ca  =  Qg.fb/a)  has  been  plotted.  A  normalizing  factor  C  has  been  added 

so  that  0  .  b2  .  IT  Q^b/a  =  l)  =  (^(b/a  *  l),  or  since  in  this  case  b2  =  a 

C  becomes  l/a2,  which  is  constant  per  definition.  Ca  increases  to  a  peak 
value  two  to  three  times  the  Ca  for  a  pure  carbon  particle,  and  then  dim¬ 
inishes  again.  This  means  that  the  water  skin  acts  in  the  direction  to 
increase  the  absorption  effect  of  the  carbon  particle.  The  absolute  values 
decrease  considerably  if  the  size  of  the  carbon  particle  itself  gets  smaller. 
Computations  ug  to  specific  ratios  b/a  of  150  show  that  the  absorption  cross 
section  Ca  =  b^  7T  *  Qa  remains  finite  and  its  obsolete  value  stays  constant 
around  the  value  for  Ca(B/a  =  20). 

An  investigation  was  made  of  the  wavelength  dependence  of  the  absorp¬ 
tion  coefficient  of  carbon  particles.  From  Fig.  1  it  can  be  seen  that  for 
4  1*5  and  for  3  <  <£-£20  the  absorption  for  carbon  particles  of  constant 
size  is  wavelength  dependent,  and  that  forOL.  L  1.5  the  absorption  is  the 
smaller,  the  larger  the  wavelength.  The  opposite  is  the  case  for  ZCCCc.  20. 

The  range  4  1.5  corresponds  for  visible  light  to  very  small  particles 
(r  <  0.1>^).  Over  the  wavelength  range  0.3  4  A <0.7,  Qq  can  vary  by  a  fac¬ 
tor  2.  This  results  in  the  phenomena  that  industrial  smoke  quite  often 
appears  red-brown  in  the  transmitted  light  if  the  particles  are  uniform 
enough  in  size  and  very  small.  As  the  particles  become  less  absorbing 
(Figs.  2  through  7),  the  slope  of  becomes  flatter  and  the  wavelength 
dependence  gradually  disappears. 
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Fig.  8.  Absorption  coefficient  for  constant  total  particle 
size  as  function  of  water  shell  thickness. 
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Fig.  9.  Total  absorption  cross  section  for 

const,  b/a  and  increasing  particle  size. 


Absorption  cross  section  for  constant  carbon 
nucleus  size  as  function  of  water  shell  thickness 


The  influence  of  vater  absorption  in  the  infrared  region  on  the  total 
absorption  of  compound  carbon- vr.t'j.  particles  must  al^o  be  studied.  In 
Table  1  the  computed  are  listed  for  a  t:tal  particle  size  b  -  6 .365^- 
(to  be  comparable  with  Shifrin's  computationt^ )  and  for  carbon- particle  radii 
between  a  s  0.171  and  0.48  for  wavelengths  from  1.35  to  10  microns.  The 
absorption  data  for  water  in  these  computations  have  been. taken  from  Dietrich® 
and  Shifrin.?  It  can  be  seen  that  the  contribution  of  the  carbon  particle  to 
the  total  decreases  with  increasing  specific  ration  b/a,  and  if  b/a  is 
larger  than  about  10,  only  about  10  percent  of  the  total  absorption  will  be 
due  to  the  carbon  nucleus.  In  other  words,  with  respect  to  absorption  proper¬ 
ties,  a  compound  carbon- water  particle  with  large  b/a  can  be  treated  in  the 
infrared  like  a  pure  water  droplet  without  introducing  a  large  error. 


Table  1. 

Absorption  coefficients  for  water  droplets  (r  = 

6.255  »)  with 

carbon  particles  of  a  = 

0.171,  0.215, 

0.368  and  0.480pu  radius. 

and  for  pure  water  droplets.* 

A.  (m) 

^a,  water  and  carbon 

Qa,  pure  water 

b/a  =  16.6 

b/a  =  29.1 

b/a  =  17 

b/a  =  13.5  b/a  =  00 

a  =  0.171 

a  =  0.215 

a  =  0.368 

a  =  0.480 

1.35 

0.0051 

0.0062 

0.0121 

0.0171 

1.5 

0.0427 

0.0440 

0.0493 

O.0537 

2.0 

0.0790 

0.0798 

0.0851 

0.0899 

3.0 

1.1846 

1.1846 

1.1846 

1.1846 

0.894 

3.2 

1.1910 

1.1910 

1.1910 

1.1910 

3.4 

1.0143 

1.0144 

1.0153 

1.0162 

0.925 

4.5 

0.5004 

0.5005 

0.5022 

0.5045 

0.500 

5.47 

0.3716 

0.3719 

0.3734 

0.3754 

6.0 

1.1311 

1.1311 

1.1312 

1.1313 

1.127 

7.0 

0.7374 

O.7376 

0.7384 

0.7395 

0.704 

10.0 

0.5751 

0.5753 

0.5757 

0.5763 

0.575 

In  Figs.  12  througn  14  the  scattering  functions  have  been  plotted  for  a 
pure  carbon  particle  for  oc  =0.6  and  oC  =  1  and  also  aL  =  50*  As  a  compari¬ 
son,  the  scattering  function  for  a  pure  water  droplet  of  the  same  size  has 
been  added.  It  can  be  seen  that  the  angular  distribution  of  the  scattered 
light  for  small  oc  -values  depends  only  slightly  on  the  refractive  index. 
However,  the  total  amount  of  scattered  light  is  larger  for  carbon  than  for 
water  particles.  The  scattering  function  for  large- size  parameter  values 
is  particularly  interesting  since,  until  recently,  no  exact  computations  for 
large  c* -values  were  available.  With  increasing  c< ,  the  convergence  of  the 
infinite  series  of  the  an  and  bn  becomes  very  poor,  and  more  and  more  terms 
have  to  be  computed.  In  addition,  the  number  of  oscillations  in  the  i^ 
and  ±2  increases  proportionally  to  oc .  As  Van  DeHulst^  has  shown,  the 

♦Values  for  water  droplets  (r  =  6.265^&)  taken  from  Shifrin. 
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Scattering  function  for  carbon  particle  with  size  parameter  Ot  -  0.6  and  X 


Fig.  l4.  Scattering  function  for  water  droplet  and 

carbon  particle  with  size  parameter  CK  -  50 
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distance  between  two  consecutive  maxima  becomes  approximately  l80°/ c*.  This 
has  been  proved  by  computations  of  Giese10  for  water  drops  (  of  =  25,  50,  100). 
It  is  also  interesting  to  notice  that  these  oscillations  for  absorbing  parti¬ 
cles  become  completely  suppressed  for  angles  9  ^  50°.  The  broad  maximum 
between  l40°  and  145°  for  water  droplets  which  indicates  the  fogbow  is  not 
present  for  absorbing  particles.  It  is  also  important  to  notice  that  for 
such  large  particles  most  of  the  total  scattered  light  is  scattered  into  the 
narrow  forward  direction  of  about  zero  to  2°,  the  intensities  per  unit  solid 
angle  being^tc  10  3  times  the  ones  for  back- scattering  angles.  Ikis  property  is 
independent  of  the  refractive  index  and  the  same  for  all  above  25.  The 
physical  explanation  for  this  is  that  for  large  particles  the  diffraction  of 
light  occurs  only  in  these  very  forward  angles.  For  most  measurements  it 
cannot  be  distinguished  from  the  direct  light,  and  therefore  the  Qs  for  such 
particles  can  be  as'sumed  as  1  instead  of  2. 

Absorption  of  Badlatlon  by  Carbon  Particles  Suspended  in  the  Atmosphere 

Before  studying  the  absorption  properties  of  a  carbon-particle  cloud, 
the  optimum  combination  of  particle  number  and  size  with  respect  to  absorp¬ 
tion  for  carbon  particles  must  be  found  for  a  given  total  carbon  mass. 

Ike  absorption  coefficient  of  a  particle  cloud  per  unit  volume  is 

</a  =  N  .  r27T«a  (°0>  (8) 

where  N  =  number  of  particles  per  unit  volume.  If  the  volume  of  one  parti¬ 
cle  is  Vp  =  47T  r3/3  and  the  total  mass  M  of  particles  per  unit  volume_is. 

Vt  •  3  (  $  =  density),  the  number  of  particles  per  unit  volume  is 
F  «  Vt  .  3/  4  TTr3  -  3^/4  T  r3  .  S*  And  ^  becomes 

^  .=  3  *L  Jb  Qa  (at),  (9) 

d  a  4  $  r 

or  with  r  =  oc  .  ^  /2  7T , 

/£.=  37rJL  - 1 -  Qa  (c<  ) 

0  a  f  2  ^  *  ot 

^  =  const  Qa  )•  (10) 

In  Pig.  15,  the  function  Qa  (c*)/°c  is  plotted,  and  the  graphical  solution 
Sives  fiaax  for  =  0<6* 

Por  sunlight  with  an  intensity  peak  at  /\  =  0.5  micron,  oC  =  0.6  corres¬ 
ponds  to  a  =  0.05  micron. 

It  should  be  mentioned  here  that  carbon  black  of  this  size  is  commer¬ 
cially  available;  however,  it  is  very  likely  that  these  small  particles  will 
form  larger  .conglomerates  during  the  dispersion  process.  Laboratory  experi¬ 
ments  on  this  problem  wiil  be  necessary. 
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.‘After  all  necessary  properties  of  single  particles  were  investigated, 
a  study  of  the  processes  of  such  particles  in  a  cloud  was  begun.  The  follow- 
ing Considerations  and  computations  are  based  mainly  on  a  thesis  by  G.  Korb 
(l9pl.  University  of  Munich,  Germany)  whose  studies  were  carried  out  on 
Signal  Corps  Contract  EUC-1612.  The  object  of  this  study  was  to  calculate 
the  absorption,  transmission,  and  reflection  of  radiation  propagating  through 
a  cloud.  In  order  to  understand  the  significance  of  the  results,  a  brief 
summary  is  given  of  Kerb's  solution. 


In  order  to  include  multiple  scattering  effects  in  the  cloud,  a  solution 
must  be  found  for  the  general  radiation  transfer  equation  (Chandrasekhar-1-'1-). 
Korb ' s  approach  is  to  solve  the  equation  for  2  +  27 T  radiation  fluxes — one 
upward,  one  downward  and  2  7T  fluxes  in  the  horizontal.  This  approximation 
reduces  Chandrasekhar's  equation  to  a  system  of  linear  differential  equations 
whose  solutions  give  the  absorption,  transmission,  and  reflection  of  radia¬ 
tion.  The  absorption  is  given  as  the  amount  of  energy  absorbed  inside  the 
cloud  by  particles  and  water  vapor.  The  method  also  takes  into  consideration 
the  absorption  of  radiation  by  water  vapor  above  and  below  the  cloud.  These 
computations  of  the  absorption,  transmission,  and  reflection  have  to  be  car¬ 
ried  out  separately  for  narrow  spectral  bands  of  the  sun  spectrum  (depending 
on  the  spectral  dependence  of  the  particle  parameters),  and  eventually  have 
to  be  added  up  to  give  the  total  values.  If  the  albedo. of  the  ground  surface 
is  0,  the  same  procedure  must  be  applied  to  so-called  apeiliary  fluxes,  i.e., 
the  diffuse  radiation  fluxes  reflected  from  the  ground  surface ;  and  the  result¬ 
ing  absorption,  transmission,  and  reflection  of  the  primary  and  secondary 
fluxes  must  be  added  up. 

* 

The  first  part  of  Kerb's  study  included  the  computation  of  absorption, 
transmission,  and  reflection  for  three  model  clouds  of  pure  water  droplets — 
a  low  cloud,  a  high  cloud,  and  a  cloud  with  large  vertical  extent. 

Korb ' s  method  was  applied  here  to  calculate  the  absorption  properties 
of  a  pure  carbon  particle  cloud,  a  cloud  mixed  of  droplets  and  carbon  parti¬ 
cles,  and  a  cloud  of  compound  concentric  carbon-water  droplets,  all  three 
clouds  having  a  vertical  extent  of  1,000  meters,  with  their  base  near  the 
ground  surface: 


(a)  Pure  carbon  particle  cloud.  The  carbon  particle  concentration  is 

lcA  particles  per  cm^,  0.1- micron-diameter  size.  This  is  equivalent  to  approxi 
mately  20  pounds  of  carbon  black  per  1  km^  for  a  1,000-meter- thick  cloud.  The 
water- vapor  content  in  the  cloud  is  assumed  to  be. 2  .  107°  g/cm3;  the  total 
water- yapor  content  above  the  cloud  (above  850  mb)  is  1.5  cm  ppw  (according 
to  an  average  distribution  over  middle  latitudes  for  summer,  Korb,  et  al^)« 


(b)  Mixed  cloud.  The  cloud  consists  of  droplets  of  4  microns  radius, 

597/cm^,  and  1CH  carbon  particles/cm3  of  0.05  micron  radius.  The  water- 
vapor  content  above  the  cloud  is  the  same  as  before--!. 5  cm  ppw;  the  satura¬ 
tion  water  vapor  density  in  the  cloud  is  assumed  to  be  8  .  10”^  g/cm3  at  about 
+7.5°C,  and  the  albedo  of  the  ground  is  zero. 


(c)  Compound  particle  cloud.  This  is  a  cloud  of  particles  which  con¬ 
sists  of  a  spherical  carbon  nucleus  with  a  concentric  water  shell.  The 
total  particle  size  again  is  four  microns  radius  and  the  concentration  is 
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597/c“3,  (This  number  results  from  a  water  content  in  the  cloud  of 
0 .26  g/m5.)  All  the  carbon  particles  of  the  cloud  in  paragraph  (a)  will 
be  inside  the  droplets,  which  gives  10,000/600  =16.7  carbon  particles/ 
droplet.  The  size  of  the  combined  carbon- particle  nucleus  becomes  approxi¬ 
mately  0.125  micron  radius,  which  gives  a  specific  ratio  b/a  =  32.  The 
water- vapor  content  in  and  above  the  cloud  is  the  same  as  in  paragraph  (b). 


The  sun  spectrum  from  0.2  through  10  microns  has  been  divided  into 
17  bands,  and  for  each  band  the  following  quantities  were  computed: 

k  _  Ka  =  £  n7Tr2  Q*  +  £  v  (A) 

Ke  n7Tr2  ^  +  ?w  *w  (A) 

and  T  =  j'  ICg  (A)  *  dh,  the  optical  density  of  the  cloud. 

t 

Then  the  intensity  of  sun  radiation  incident  on  the  cloud,  EQ  (  A ) 

was  computed;  and  finally  the  absolute  values  of  absorbed,  transmitted, 
and  reflected  energy  were  obtained.  All  computations  have  been  carried  out 
for  zenith  distances  of  the  sun  of  0°,  30®,  60°  and  80°. 


(ID 

(12) 


The  results  of  the  computations  are  plotted  in  Fig.  16.  In  the  graphs 
for  a  mixed  cloud,  the  absorption  curve  for  a  pure  water-droplet  cloud 
(taken  from  Korb's  report)  has  been  added  for  comparison.  Also  added  is  the 
curve  for  Eq,  i.e.,  the  intensity  of  energy  incident  on  the  cloud  top  after 
passing  through  the  upper  atmosphere.  Approximately  0.2  •  cal/ cm2  •  “in- 
are  absorbed  by  water  vapor  in  the  upper  atmosphere  above  the  cloud,  for 
vertical  incidence.  In  the  case  of  the  pure  carbon- particle  cloud,  the  per¬ 
centage  of  absorbed  and  transmitted  energy  is  almost  equal,  with  very  little 
energy  (<  5  percent)  reflected  back  from  the  cloud  top.  FOr  very-low  sun 
elevations,  such  a  cloud  absorbs  around  80  to  90  percent  of  the  incident 
radiation,  although  the  absolute  amount  of  absorbed  energy  becomes  very  small 
for  low  sun  elevations.  If,  as  a  comparison,  the  absorption  in  such  a  pure 
carbon  cloud  for  single  scattering  only  is  calculated  applying  Bouguer's  law, 


Ea 


(1  -  e' 


-Ja«L 


), 


(13) 


substitution  is  made  for  =  N  •  =  5-5  '  10"7  cm"  >  and  the 

path  length  for  zenith  distance  z  =  0  is  10^  cm.  Then,  Ea/E0  =  5-5  percent. 

The  absorption  by  water  vapor  in  the  cloud  is  approximately  2.5  percent. 

This  would  give  a  total  absorption  of  about  7  percent,  only  one-sixth  of 
the  actual  absorption. 

The  optical  properties  of  the  mixed  cloud  are  considerably  different. 
About  60  percent  of  the  incident  radiation  is  reflected  from  suck  a  cloud 
so  that  only  a  small  portion  of  tne  radiation  penetrates  into  the  cloud. 
However,  it  can  be  seen  that,  because  of  the  carbon  particles,  about  three 
times  more  energy  is  absorbed  than  in  the  pure  water-droplet  cloud.  Only 
about  12  percent  of  the  incident  light  is  transmitted  through  this  cloud. 
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Fig.  16.  Computed  transmitted  energy  (Ed),  absorbed  energy  (Ea)  and  reflected 

energy  (Eg)  for  three  cloud  models,  as  function  of  sun  zenith  distance. 


A  very  interesting  result  is  that  the  optical  properties  of  the  mixed, 
and  the  compound  cloud  are  almost  equal,  the  absorption  of  the  mixed  cloud 
oeing  about  2  to  2.5  percent  (of  eI)  higher  than  for  the  compound  particle 
cloud.  From  this?  it  can  be  concluded  that  the  form  in  which  the  carbon  par¬ 
ticles  are  distributed  in  the  water-droplet  cloud— whether  they  are  in  between 
the  droplets,  or  inside,  or  on  the  surface- -makes  little  difference  in  the 
absorption  coefficient  so  long  as  the  total  carbon  mass  is  about  the  same. 

This  result  is  not  too  surprising  if  one  looks  at  the  spectral  distribution 
of  the  absorbed,  energy.  About  k-5  to  50  percent  of  the  total  absorbed  energy 
is  absorbed  in  the  infrared  between  two  and  ten  microns  wavelength,  although 
there  is  only  about  27  percent  of  the  total  incident  energy  in  this  wavelength 
range.  This  is  because  of  the  high  absorption  of  water  in  the  infrared.  This 
wavelength  region  is  therefore  influenced  very  little  by  the  carbon  absorption. 


In  the  visible  range.  Figs.  10  and  11  show  that  for  large  b/a  the  absorp¬ 
tion  cross  section  Ca  =  b^  TT  *  Qg  remains  rather  constant  and  is  only  slightly 
higher  than  the  one  for  the  carbon  nucleus  without  water  shell.  This  result 
is  important  because  it  is  not  actually  known  what  happens  to  the  carbon  par¬ 
ticles  in  the  cloud  after  seeding — whether  they  will  attach  themselves  to  the 
surface  of  droplets,  or  whether  most  of  them  will  remain  suspended  in  the  air. 
It  can  be  assumed  that  the  rather  special  case  of  concentric  compound  parti¬ 
cles  will  exist  very  seldom — probably  only  if  condensation  of  water  on  the 
carbon  particles  occurs,  most  likely  by  radiational  cooling  of  the  carbon 
particles  during  nighttime. 


So  far  in  these  calculations  it  has  been  assumed  that  the  albedo  of  the 
ground  surface  is  zero.  This  means  that  the  upward-directed  radiation  flux 
at  the  bottom  of  the  cloud  is  zero.  dhe  more  light  there  is  reflected  from 
the  ground  back  into  the  cloud,  the  higher  the  total  absorption  in  the  cloud 
will  be.  Korb  has  computed  the  absorption,  transmission,  and  reflection  of 
pure  water-droplet  clouds  for  various  albedo  values  between  zero  and  0-9,  and 
from  his  data  it  can  be  estimated  that  in  the  case  of  a  mixed  and  compound 
particle  cloud  the  absorption  would  be  increased  by  about  2  to  3  percent 
(of  Eo)  if  the  ground  albedo  were  0.9  instead  of  zero.  The  increase  in 
absorption  is  so  small  because,  in  addition  to  the  ground-reflected  radiation,  a 
great  percc  lege  of  radiation  is  reflected  down  from  the  bottom  surface  of  the 
cloua.  ZLilz  ircreases  the  trar«missivity  of  the  cloud  considerably.  For  the 
carbon  particle  cloud,  of  course,  higher  ground  albedo  would  result  in  a  much 
higher  absorption.  The  denser  a  cloud,  the  less  influence  the  ground  albedo 
has  on  the  absorption  of  the  cloud. 


For  energy  processes  inside  the  cloud,  the  vertical  distribution  of  the 
total  absorbed  energy  within  the  cloud  must  be  considered.  From  Korb's  thesis 
it  follows  that  for  a  cloud  with  optical  thickness  30  km"-*-  (as  in  the  compound 
and  mixed- cloud  model),  for  zenith  distance  z  =  30°  and  ground  albedo  As  =  0 
and  Kg/Kg  =  0.003,  50  percent  of  the  total  absorbed  radiation  is  absorbed  in 
the  upper  300  meters  of  the  cloud;  whereby,  for  2  4  T430,  the  intensity  of 
the  downward- directed  radiation  flux  almost  follows  a  distribution  of  the 
form  E  =  -a  'L  +  b.  The  denser  the  cloud,  the  greater  the  percentage  of 
absorbed  energy  in  the  upper  portion  of  the  cloud. 


If  results  of  this  study  are  compared  with  previous  calculations  by 
other  investigators,  the  rather  detailed  study  by  Smith,  et  al  ^  should  be 
mentioned.  They  have  computed  the  transmittance,  reflectance,  and  absorption 
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of  a  mixed  cloud  of  droplets  and  small  absorbing  particles,  neglecting  the 
absorption  in  droplets  and  applying  the  two-flux  model  by  S.  Fritz. For 
high  optical  thicknesses  the  two- flux  model  gives  about  10  percent  lower 
transmission  values,  according  to  comparisons  in  Korb's  report.  Using 
calculations  of  Smith,  Wexler  and  Glasser,  ^  a  cloud  transmission  of  about 
5  to  10  percent,  and  a  reflectance  of  about  80  percent,  a  10  to  15  percent 
absorption  is  obtained  for  the  compound-particle  cloud.  This  is  for  verti¬ 
cal  incidence.  Computations  made  under  this  study  give  about  25  percent 
absorption,  60  percent  reflection,  and  15  percent  transmission.  It  seems 
that  the  two- flux  model  overemphasizes  the  reflectance  of  a  cloud.  For 
detailed  discussions  of  the  accuracy  and  significance  of  data  obtained  with 
the  2  +  2 Iff lux  model,  reference  is  made  to  Korb's  thesis. 

Finally,  an  estimate  is  made  of  how  much  cloud  evaporation  or  dissipa¬ 
tion  can  be  achieved  with  this  absorbed  energy.  The  radiation  energy  which 
is  absorbed  within  the  cloud  is  used  for  heating  of  the  total  air  mass  and 
water-droplet  mass,  and  then  for  evaporation  of  the  droplets.  If  it  is 
assumed  that  the  absorption  and  evaporation  process  occurs  slowly  enough  so 
that  equilibrium  conditions  are  close  between  all  the  components  and  the 
atmosphere,  then,  for  the  energy  budget, 

dq  -  cpa  •  S'a  *  dT  +  cy  •  w  •  '  dT  +  L  •  w.  (it) 


cpa  =  specific  heat  of  air  for  const,  pressure,  0.2k  cal/g  °C 
3  a  =  density  of  air,  at  7-5°  and  90°  mb  =  1.12  g/cm3 
Cy,  =  spec,  heat  and  density  of  water; 
w  =  liquid  water  content  of  cloud  g/m°, 

L  =  latent  heat  of  vaporization  of  water,  590  cal/g  °C. 

Substituting  these  values  into  equation  (lk), 

dq  =  269  +  0.155  +  87  (cal/m3) 

=  356  cal/m3  for  a  1°C  temperature  increase. 


Heating  of  the  carbon  particles  has  been  neglected  completely  because  their 
mass  is  only  about  l/lOO, 000th  of  the  water  content  and  it  is  seen  that  even 
the  energy  required  for  the  heating  of  the  droplets  is  negligible.  The  tem¬ 
perature  rise  which  is  required  to  increase  the  saturation  water  vapor  con¬ 
tent  from  the  original  8  g/m3  to  8.15  g/m3  is  approximately  0.3°C.  From 
equation  (lk),  therefore,  is  obtained  the  energy  required  to  evaporate  the 
compound-particle  cloud  and  to  achieve  the  necessary  temperature  increase  of 
0.3°C  to  about  100  cal/m^. 

In  the  case  of  the  mixed  cloud  and  the  compound-particle  cloud,  the 
absorbed  energy  for  sun  zenith  distances  between  zero  and  30°  is  about 
O.k  cal/cm^  min.  The  air  column  below  1  cm^  is  1  •  105  cm3  =  0.1  m3,  which 
gives  an  average  absorbed  energy  of  k  cal/m3  •  min.  This  vould  result  in  a 
necessary  absorption  time  of  about  25  minutes.  For  a  sun  zenith  distance  of 
60 °,  the  absorbed  energy  would  be  only  about  0.2  cal/cm^  min.  and  the  absorp¬ 
tion  time  therefore  twice  as  long--apprcximately  50  minutes. 
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Actually,  it  must  be  expected  that  the  upper  100  meters  of  the  cloud, 
where  the  absorption  per  unit  volume  is  about  twice  the  average  absorption, 
will  start  to  warm  up  and  evaporate  first.  This  will,  again,  allow  more 
radiation  to  penetrate  into  the  lower  layers  of  the  cloud.  However,  as  soon 
as  the  upper  layer  starts  to  warm  up  it  will  also  begin  to  rise  and  cool  off 
again.  It  certainly  will  be  necessary  to  provide  more  than  the  minimum  energy 
for  a  0.3°C  temperature  increase  if  the  evaporation  process  is  to  be  maintained 
long  enough  to  evaporate  the  cloud  completely. 

Van  Straten^  suggested  as  a  hypothetical  process  for  the  cloud  dissipa¬ 
tion  an  evaporation  and  condensation  process  between  droplets.  Droplets  con¬ 
taining  many  carbon  particles  would  heat  up  more  than  pure  water  droplets, 
and  therefore  moisture  would  evaporate  from  the  warmer  droplets  with  the 
carbon  particles  onto  the  colder  pure  water  droplets. 


A  droplet  containing  carbon  particles  will  absorb  energy  and  conse¬ 
quently  its  temperature  will  rise.  This  again  will  cause  an  energy  loss  by 
conduction  and  evaporation.  Hie  energy  budget  can  be  written  in  the  follow¬ 
ing  way: 


<iEa  =  c  •  m  •  dT  +  r  •  K  •  (T  -  T0)  +  ^Tr2  •  f  •  L  -4£_  ^ 

at"  at  at 

c  =  specific  heat  of  water  =  1  cal.  degr.-^-  *  g"d 
m  =  mass  of  droplet 

r  =  droplet  radius  _ 

K  =  thermal  conductivity  of  air  =5*8*  lO-5  cal.  cm 
sec- a  •  degr'l 
T  =  droplet  temperature 
T0=  air  temperature 
$  =  water  density  =  1  g  .  cnf 3 

L  =  latent  heat  of  evaporization  =  590  cal  .  g”1  .  degr"1 

whereby  ar  =  D -  ($d  -£y)  (16) 

dt  J  '  r 

D  =  coefficient  of  diffusion  of  water  vapor  =  0.225  cm^  .  sec~d 

Jd,  =  vapor  density  of  water  at  the  droplet  surface  and 
far  away  from  the  droplet,  respectively. 


§a  is  a  function  of  T  and  r,  according  to  Kelvin's  formula: 


■? 


d 


.  pw(T)  •  exp( 


M  | 
R’r’T 


(17) 


6~~  =  specific  surface  energy  of  water  =  7^.5  cal  •  cnf^ 
M  =  molecular  weight  of  water  =  18 
R  =  universal  gas  constant  =  I.98  cal  •  g"d  •  degrT^. 


Substituting  (16)  and  (17)  in-  (15), 


dE 

_ 1  =  c  * 

dt 


dT  +  hTTr  •  k(T-T0)  +  U^rLD 
dt 


(yw  (T)J  •  exp(^|!)  =  ^w(To).  (18)' 


m 
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As  a  first  approximation,  it  is  assumed  that  r  is  constant  and  the 
following  abbreviations  are  used: 

Of  =  c  •  m  =  2.68  •  lO”10  cal  •  degr”"* 

^  =  4  7T  rK  =  2.96  •  10-7  caj_  .  degr”1  •  sec  1 

^"=4  rLD  =  0.67  cal  g  cm"3  sec"1 

c  -  =  0.082  degr. 

0  R  •  r 


dEa 

IT 


q.  = 


r2  Tf  ’  Qa  *  Eo 
60 


215  ♦  10“ 11  cal 


sec" 


1 


with  Qa  =  0.002,  Eq  =  1«5  cal 


cm 


-2 


min 


-1 


We  then  get  g 

q  =  Q|  +  0(T-To)  +  (T)  eT  -Sy  (TQ)j  , 

£ 

°r  *§+  ot +a*  (t)  *T = q +(3to + i^w  (tc) 


(20) 


c  =  q  +  /3tq  +  (f  y  (T  )  is  constant.  The  magnitude  of  each  of  the  three 
quantities  (a  =  2.5  •  1CT11,  (STq’  =  8*3*  10"5,  and  ^  •  S*W(T0)  = 

5-4  •  10"°),  shows  that  q  will  have  to  be  100,000  times  bigger  to  be  of  any 
influence  on  the  result  of  this  equation.  This  confirms  Smith,  Wexler,  and 
Glaser's  conclusion12  that  lCr  particles  per  droplet  would  be  necessary  to 
achieve  a  sufficient  rise  of  temperature  in  the  droplet. 


An  analytical  integration  of  equation  (18)  will  be  extremely  difficult, 
if  not  impossible.  However,  since  in  this  study  there  is  not  much  interest 
in  the  transient  phase,  dT/dt  =  0  is  assumed  and  T  is  calculated  only  after 
reaching  equilibrium  conditions. 

£ 

Q  T  +  fSv  (T)  eT  -  c  =  0.  (21) 

The  numerical  solution  gives  T  =  2Q0.h6°K  or  7-3°C.  The  temperature  of  the 
particle  has  decreased  by  0.2°C  from  the  initial  temperature,  T0  =  W 
From  an  assumed  air  temperature  of  7.5°C,  evaporation  .’ue  to  the  curvature 
of  the  droplet  was  obtained.  Also,  from  equation  (l8),  it  can  be  seen  that 
the  condensation  term  is  of  the  order  10"7  and  the  evaporation  term  is  of 
the  order  10"°;  and  the  first  term,  the  internal  energy  increase,  is  of  an 
order  <  10"1(^,  which  means  that  most  of  the  absorbed  energy  is  lost  by  con¬ 
duction  to  the  surrounding  air.  This  loss  of  energy  to  the  air  will  actually 
be  even  greater  since  the  heat  transfer  does  not  occur  only  by  conduction 
but  also  by  convection.  Consequently,  the  air  temperature  will  rise  almost 
as  fast  as  the  particle  temperature.  It  has  been  shown  that  for  small 
absorption  rates  the  particle  temperature  may  be  even  below  the  air  tempera¬ 
ture  because  of  evaporation  due  to  vapor- pressure  difference.  Most  of  the 
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evaporated  water  from  the  droplets  will  be  needed  to  maintain  saturation  in 
the  heated  air,  and  only  a  very  small  portion  would  be  available  for  conden¬ 
sation  on  pure  droplets.  Therefore,  it  cau  le  stated  definitely  that  an 
evaporation- condensation  process  cannot  be  conf!  dered  as  a  possible  cloud- 
dissipation  process,  and  therefore  the  only  possible  dissipation  process  for 
carbon- seeded  clouds  will  be  evaporation  due  to  heating  of  the  total  air  mass 
which,  as  has  been  shown,  requires  very  large  quantities  of  carbon  black. 

This  report  has  been  limited  to  absorption  processes  only,  and  no  dis¬ 
cussion  has  been  made  of  any  emission  effects  from  carbon  seeded  clouds. 

For  more  information  on  the  energy  budget  in  clouds,  a  study  is  certainly 
necessary.  As  a  result  of  the  very  encouraging  results  of  the  first  part 
of  Korb's  study  and  the  results  reported  here,  the  studies  by  the  University 
of  Munich,  Germany,  have  been  extended  under  Signal  Corps  contract  to  also 
cover  various  models  of  pure  and  mixed  carbon  and  water  clouds  to  a  much 
hx0ur"‘  extent  and  accuracy  than  has  been  done  in  work  reported  here. 

This  extender  3tudy  will  also  include  energy  losses  due  to  emission  from  the 
cloud. 

CONCLUSIONS 

A  summary  of  the  results  indicates  that  for  a  model  cloud  of  approxi¬ 
mately  600  droplets/ cm^,  4  micron  radius,  from  ground  to  1,000  meters  alti¬ 
tude  (stratus- type  cloud),  seeding  with  20  pounds  of  0.05  micron  radius 
carbon  particles  will  increase  the  absorption  of  sun  radiation  within  the 
Cloud  from  about  0.15  to  0.4  cal/cm2  •  min  for  0°  to  30°  sun  zenith  distance. 

The  absorption,  transmission  and  reflection  are  almost  equal  for  a 
mixed  cloud  of  water  droplets  and  carbon  particles  and  a  cloud  of  compound 
concentric  carbon- water  particles.  This  indicates  very  strongly  that  for 
absorption  studies  it  is  not  very  critical  to  know  in  what  form  the  carbon 
particles  are  distributed  within  the  cloud,  i.e.,  whether  they  are  contained 
vithin  the  cloud  droplet  or  attached  to  its  surface. 

From  estimates  of  the  energy  required  for  evaporation  of  the  model 
cloud  it  is  concluded  that  for  realistic  field  experiments  it  will  be  nec- 
tssary  to  use  seeding  rates  of  the  order  of  at  least  100  pounds  of  carbon 
black  per  km2  during  high  sun  elevations  and  for  relatively  stable  cloud 
formations.  The  necessary  seeding  rate  may  go  up  to  several  hundred  pounds 
per  km2  under  less  favorable  conditions. 

An  analysis  of  the  energy  budget  in  the  seeded  cloud  yields  that  an 
evaporation- condensation  process  will  not  be  possible,  and  that  evaporation 
Of  a  cloud  due  to  heating  of  the  total  air  mass  will  be  possible  only  under 
very  favorable  conditions. 
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